Dynamic instability of microtubules: effect of catastrophe-suppressing drugs 
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Microtubules are stiff filamentary proteins that constitute an important component of the cy- 
toskeleton of cells. These are known to exhibit a dynamic instability. A steadily growing micro- 
tubule can suddenly start depolymerizing very rapidly; this phenomenon is known as "catastrophe" . 
However, often a shrinking microtubule is "rescued" and starts polymerizing again. Here we de- 
velope a model for the polymerization-depolymerization dynamics of microtubules in the presence 
of catastrophe- suppressing drugs. Solving the dynamical equations in the steady-state, we derive 
exact analytical expressions for the length distributions of the microtubules tipped with drug-bound 
tubulin subunits as well as those of the microtubules, in the growing and shrinking phases, tipped 
with drug-free pure tubulin subunits. We also examine the stability of the steady-state solutions. 

PACS numbers: 87.10.+e, 82.35.Pq, 87.15. Rn 
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I. INTRODUCTION 

Microtubules are filamentary proteins and constitute 
a major component of the cytoskeleton of the eukary- 
otic cells [l|, Hf. The dynamic cytoskeletal scaffolding 
not only supports the cell architecture and gives rise 
to changes in the shape of the cell but the network of 
its constituent filamentary proteins also provides path- 
ways for intra-cellular transport. In other words, a wide 
range of dynamical processes, which are essential for 
sustaining life, are driven by the dynamic cytoskeleton. 
Therefore, a clear theoretical understanding of the fun- 
damental physical principles behind the polymerization- 
depolymerization dynamics of the microtubules is ex- 
pected to provide deep insight into the physics of cell 
shape transformations, cell motility, etc. as well as mech- 
anisms of many sub-cellular processes like, for example, 
chromosome segregations during mitosis (i.e., cell divi- 
sion). 

Dynamic instability [H, HJ is one of the unusual non- 
equilibrium processes that dominate the dynamics of 
microtubule polymerization. Each polymerizing mi- 
crotubule persistently grows for a prolonged duration 
and, then makes a sudden transition to a depolymeriz- 
ing phase; this phenomenon is known as "catastrophe" . 
However, the subsequent rapid shrinking of a depolymer- 
izing microtubule can get arrested when it makes a sud- 
den reverse transition, called "rescue" , to a polymerizing 
phase. The elongation of a microtubule takes place by 
a reversible, non-covalent attachment of a tubulin dimcr 
from the tubulin solution. It is now generally believed 
that the dynamic instability of a microtubule is triggered 
by the loss of its guanosine triphoshate (GTP) cap be- 
cause of the hydrolysis of GTP into guanosine diphos- 
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phate (GDP). But, the detailed mechanism, i.e., how the 
chemical process of cap loss induces mechanical instabil- 
ity, remains far from clear. 

The dynamics of polymerization-depolymerization of 
microtubules and the phenomena of "catastrophe" and 
"rescue" [!, H, 0, II], have been studied extensively 
over the last decade using simple theoretical models 

[j| M, ED, Hi HI Hi, HE HE EE m HI H El- o ne of 

the earliest models of polymerization-depolymcrization 
dynamics of microtubules was proposed by Hill |9| and 
subsequently extended by Rubin and coworkers [lO . [l9| . 
The length of the microtubules is discrete in the Hill 
model which is formulated in terms of two infinite sets 
of coupled ordinary differential equations. Dogterom and 
Leibler [HI [12], however, treated the length as a contin- 
uous variable and the dynamical equations were reduced 
to just two coupled partial differential equations. 

A large number of different types of antimitotic 
drug molecules are known to bind with free tubu- 
lins in solution and/or with tubulins in microtubules. 
The polymerization-depolymerization dynamics of mi- 
crotubules, which play crucial roles in mitosis, is strongly 
influenced by these drugs. The effects of various drug 
molecules, e.g., colchicine, paclitaxel, vinca alkaloids and 
taxol, etc., on the dynamics of microtubules have been 
investigated experimentally for several years [22l . l23l |24| , 
partly because of their potential clinical use in combating 
cancer [H, [H, These drugs can be broadly classi- 
fied into two groups. One group consists of microtubule- 
destablizing agents whereas the members of the other 
group are microtubule- stabilizing agents. 

In this paper we are interested in the effects 
of microtubule-stabilizing agents (i.e., catastrophe- 
suppressing drugs) on the length distributions of the mi- 
crotubules in the absence of any GTP and GDP in the 
system. The generic drug molecules of our interest are 
assumed to bind rapidly with free tubulins in solution; 
when such a tubulin-drug complex binds with the grow- 
ing end of a microtubule, the drug-capped microtubule 
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gets stabilized because of the strong suppression of catas- 
trophe phenomenon [26|, [27| • 

Recall that there are four main rate constants 
(or, frequencies) that characterize the four impor- 
tant processes involved in microtubule polymeriza- 
tion/depolymerization dynamics. Two of these are the 
rate constants for the attachment and detachment of 
tubulin dimers in the polymerization and depolymeriza- 
tion phases, respectively while the remaining two are 
the frequencies of catastrophe and rescue. We assume 
that the generic drug has the following two effects on the 
polymerization dynamics: (i) it reduces the frequency of 
catastrophe to such a large extent that the microtubules 
capped with drug-bound tubulins do not exhibit catas- 
trophe at all, and (ii) it also affects the rate of elongation 
of the microtubules because the rate constant for the at- 
tachment of a drug-bound tubulin is, in general, differ- 
ent from that of a drug-free tubulin. The effects of real 
catastrophe-suppressing drugs are quite complicated and 
depend also on the dosage of the drug; some even induce 
a "paused" phase [2^ |. 

The aim of this paper is to investigate theoretically 
the generic effects of catastrophe-suppressing drugs by ex- 
tending the earlier theoretical models, developed by Hill 
@ and Freed [2(| , for the dynamics of polymerization of 
drug- free pure tubulins. We derive exact analytical ex- 
pressions for the steady-state distributions of the lengths 
of the microtubules tipped with the drug-bound tubu- 
lin subunits as well as those of microtubules tipped with 
pure (i.e., drug- free) tubulin subunits. We carry out lin- 
ear stability analysis of the steady-state distributions and 
physically interpret the implications of the spectrum of 
the eigenvalues of the stability matrix. 

This paper is organized as follows. In section [TT] we 
briefly review some of the relevant earlier theoretical 
models of dynamic instability. In paricular, we sum- 
marize the mathematical frameworks of the Hill model 
Q and the recent Freed model 20] of dynamic insta- 
bility of microtubules. In section IIIII we propose an ex- 
tension of the Hill model so as to capture the effects of 
the catastrophe-suppressing drugs in a simple way. The 
model is made more realistic in section IIVI by treating 
the dynamics of the concentration of the drug-free tubu- 
lins explicitly following a Freed-like approach. The paper 
ends with a conclusion section fVl 



II. BRIEF REVIEW OF EARLIER MODELS OF 
DYNAMIC INSTABILITY 

A. Hill model 

Microtubules consist of 13 protofilaments, each con- 
sisting of monomeric units (actually a a — (3 hetero- 
dimer) each of which is approximately 8 nm long. On 
the other hand, the polymers in the Hill model are one- 
dimensional. Therefore, some authors (see, for example, 
(28|) identify "monomeric units" of the one-dimensional 



Hill model to have a length of approximately 8/13 = 0.6 
nm. 

Let Pft(n,t) and P^(n,t) be the probabilities of find- 
ing a microtubule of length n, at time t, in the growing 
(+) and shrinking (— ) phases, respectively. Moreover, let 
P#(0,i) be the probability that the MT nucleating site 
is empty at time t. Furthermore, we denote the growth 
rate of the growing microtubules and decay rate of the 
shrinking microtubules by p^ andp^, respectively, while 
the frequencies of catastrophes (the transition from grow- 
ing to the shrinking phase) and rescues (the transition 
from the shrinking to growing phase) are denoted by the 
symbols and p_,, respectively. Interestingly, the 
four parameters p^ , p^ , p+_ and p H _ + were measured 
as functions of free tubulin concentration first in 1988 
[29( by observing single microtubules using video light 
microscopy. However, the concentration dependence of 
these parameters, if any, does not appear explicitly in 
the Hill model. 

The Master equations for these probabilities are given 

by 1 



alt 



:pfP+(n-l,t)+p H + P H (n,t)- 

(pf +P+-)P&(n,t), forn>2, (1) 



dt 



p»P„(n + l,t)+p»_P+(n,t)- 
(pf +p H + )P H (n,t), forn> 1, (2) 



^^=P»PH(0,t)+p" + P H (l,t)-(p^+p»_)P+(l^ 

(3) 

and 

^m = ^ PHM+P H P - {ljt) _ (4) 

Imposing the normalization 

CO oo 

]T^H + ]TP H (n) + P ff (0) = l, (5) 

n—1 n— 1 

the steady-state solutions of the equations CD)-® are 
given by 



and 



with 



P+(n)=x n H P H (0), 
p-(n) = x n H - 1 y H P H (0) 



Ph(0) 



1 + Vh 



x H 



jff (iff 

Pd(Pg +P+-) 



(0) 
(7) 

(8) 
(9) 
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and 



Pg 

VH = -jY- 

Pd 



(10) 



In order that the distributions P^(n) are decreasing, 
rather than increasing, functions of n we must have xh < 



1 which imposes the constraint p^ (p d - 
P+-), i-e-, 



pf P?+ < Pfp%- 



<Pd(Pg~ 



(11) 



on the magnitudes of the parameters. Interestingly, in 
the Dogterom- Leibler model [III > a steady-state charac- 
terized by exponentially decaying distributions P^ of the 
lengths of the microtubules is attained only if the con- 
dition (fTTj) is satisfied by the parameters; otherwise, the 
system never reaches any steady-state and mean of the 
(Gaussian) distribution of the lengths of the microtubules 
continues to increase linearly with time. 



B. Freed model 

It has been realized for quite some time [ljl [l?], that 
the rate of growth of the growing microtubules should de- 
pend on the availability of free tubulin monomers in the 
solution. However, in the Hill model 0] the kinetic rate 
equations do not involve the concentration of the tubulin 
monomers. Recently, Freed [2(| has generalized the Hill 
model to incorporate the dependence of the rates on the 
tubulin monomer concentration. Hammele and Zimmer- 
mann [21| carried out independent analytical calculations 
of the same phenomenon by extending the Dogterom- 
Lcibler model. Since our calculations in the section 
IIVI are based on an extension of Freed's model, we sum- 
marize here the main points of this approach. 

Suppose, po is the initial concentration of the tubulin 
subunits and p is the corresponding instantaneous con- 
centration at time t. Similarly, N$ and N are the initial 
and instantaneous concentrations of the free (i.e., with- 
out bound tubulin) nucleating sites, respectively. The 
symbols P F {n,t) in this section denote the concentra- 
tions, rather than probabilities, of microtubules in the 
growing and the shrinking phases, respectively. More- 
over, binding of a tubulin subunit with a free nucleating 
site takes place at a rate p F . Using these quantities, the 
kinetic rate equations in the Freed model [20l | can now 
be written as 



dP+{n,t) _ , 



dt 



= pi P P+ {n - 1, t) + p%P F (n, t) - 

{p F gP + plJ)P${n,t), forn>2, (12) 



dP F (n,t) 
dt 



p F P F {n+l,t)+p F _P+{n,t)- 

(p F +P F +)P F (n,t), forn>l, (13) 



dt 



ripN+pt+Pjil, t) - (p^p+p^_)P+(l,t), 

(14) 



and 



dp 



f t = -p F pN-p F pJ2PMn,t)+p F J2 P F(^)- (15) 



Moreover, tubulin mass conservation imposes the condi- 
tion 



where 



Furthermore, 



where 



Po = P + Qf + Qf 



Q F = J2nPp(n,t). 

n=l 



N Q = N + P++P F 



(16) 



(17) 



(18) 



(19) 



There is one-to-one correspondence between the pa- 
rameters and dynamical variables in the Freed model and 
those in the Hill model. For example, p F p,p F ,p F - and 
p F _ + correspond to p^ ,Pj ,p+_ and p F [ +1 respectively. 

The steady-state solution of this system of kinetic 
equations is given by [3(| 



P+ =p F N p(p F + +p F )/D 



(20) 



and 



P F =p F N p(p F _+ P F p)/D (21) 

P F (1) = f3'=p F N p(p F p F _ -p F g pp F + )/{p F d D) (22) 
where 

D = PnP F p 2 + (PdP F +PnP+--p F p F +)P+(PdP I t-+P F pP-+) 

(23J 

Finally, using a generating function technique, Freed [201] 
derived the analytical expressions 

P+(n) = ( a V) ( "- 1)/2 [(/' + /^')t4-i(A') 



-(a'cT 1/2 a'fU n - 2 (\')} (24) 



and 



P F (n) = {a'd)W 2 \pU n -i{>!) 
+ (a'c')-^ 2 (b'f -c7?')[4- 2 (A')] for n > 2, (25) 

where 

a' = (p F +p F + )/p F , (26) 
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b' = -vl-lvl, (27) 

c' =p F p/{p F p + p F + _), (28) 

d' = P F + /(p F P + P F -), (29) 

f =p F pN/{p F p + P F + _). (30) 

A' = (a' + b'd 1 + c')/[2(a'c') 1/2 ] (31) 

and U n {\') are the Chebyshev polynomial of the second 
kind given by 

U n (\') = sin[(n + 1) arccos A']/(± - A ,2 ) 1/2 , (32) 

together with t/-i(A') = 0. 

The distributions Pp{n) will be decreasing functions 
of n provided the condition a'cJ < 1 is satisfied; this 
condition imposes the constraint 

F F F F /on\ 

PgPP-+ <P d P+- ( 33 ) 

on the magnitudes of the parameters and the value of 
the drug-free tubulin subunits p in the steady-state. The 
condition (|33|) becomes identical to (fTTj) if we identify 
p F p with pf. 

While expressing the steady-state distributions Pp (n) 
in terms of Chebyshev polynomial, Freed [2(| implicitly 



assumed that |A'| < 1. However, as shown in appendix^ 
A' is, in general, larger than unity. We revise the Freed's 
result in section Hvl bv taking A' > 1 and, consequently, 
we get a different polynomial instead of the Chebyshev 
polynomial given in (|32p . 

Freed [2(| derived the exact form of the stability ma- 
trix, we define 

A± = £*P±(n) (34) 
n=l 

and the column vectors 




The, the equations obtained from the linear stability 
analysis above can be written as 

= M P V(t) + N F SP F (1) (37) 

at 

where the matrix Mp is given by 



M F 




P-+-Pn[p]ss M i3 

-p-+ o 

Pn[p}ss+Pd M33 



(38) 



with 

M13 = P n{N - [P+}ss - [P-}ss} (39) 

and 

M 33 = (p» - Pg )[P+}ss ~ Pn{N - [P_]„) (40) 

The nature of the stability of the steady-state is deter- 
mined by the eigenvalues of the matrix Mp. The steady 
state is stable if all the eigenvalues are real and nega- 
tive. On the other hand, if some roots are positive, these 
would indicate unbounded growth and the correspond- 
ing steady-state would be unstable. If the characteristic 
equation has a pair of complex conjugate roots then the 
system will either oscillate about the steady state (if the 
real part is negative) or exhibit oscillatory unbounded 
growth (if the real part is positive). 



We have carried out numerical analysis of this stability 
matrix and obtained the eigenvalues for several sets 
of values of the model parameters; the eigenvalues 
corresponding to five different values of p, for a fixed 
set of values of the other model parameters, are shown 
in the table fl] Thus, the steady-state distribution is 
stable over the entire range of p for the chosen set of 
parameters values. 



III. HILL-LIKE MODEL WITH DRUGS 

Let Pc be the rate of growth of a microtubule by 
addition of a drug-bound tubulin. Since addition of 
catastrophe-suppressing drugs to the system strongly re- 
duce the catastrophe frequency, we assume that the drug 
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p 


mi 


777,2 


777-3 


0.4 


-49.984 


-0.282 


-23.719 


0.8 


-99.999 


-0.295 


-22.507 


1.6 


-199.999 


-0.299 


-20.456 


2.0 


-250.000 


-0.299 


-19.566 


2.5 


-312.500 


-0.300 


-18.577 



TABLE I: The eigenvalues of the linear stability matrix Mf 
in the Freed model. The other common parameters for above 
table are p a = 125, p d = 900, p-+ = 0.08, p+- = 0.22, An = 
0.2 (in respective units). 



lis such that it arrests catastrophe. In other words, a 
microtubule tipped with a drug-bound tubulin can grow 
but cannot shrink. We shall use the symbol Tlh(n,t) to 
denote the probability of a microtubule, tipped with a 
drug-bound tubulin, that has length n at time t. As a 
consequence of our assumption, we do not need to con- 
sider the two quantities n^(n, t) and n^(n, t) separately; 
IL^(n,t) = for all n at all t. However, even in the 
presence of such drugs in the system, catastrophes can 
take place in microtubules tipped with drug-free tubu- 
lins. The distributions of the microtubules tipped with 
drug-free tubulin subunits in the growing and shrinking 
phases are denoted by P^(n,t) P^(n,t), respectively. 



P n 



Pnc 



o 



IOCO® 

> 



5 o 

p © 



IOOCXD 



P-+ 

IOOOO : 



I': 



-O 



I M TOC ( N : Drug-free tubuline £^ Drug-bound tubuline 
( ) Tubuline constituting MT but not located at the tip 



FIG. 1: A schematic description of our Hill-like model in the 
presence of catastrophe-suppressing drugs. The superscript 
h has been dropped from the symbols denoting the model 
parameters. 



( "' /0 = P h d P t ;(n + l,t)+p h + _P+(n,t) 



dt 
dt 

dP h (0,t) 



-(p h d +p^+)P f ;(n,t), for n > 1, (42) 



^P, l (0,t)+^ + P ft -(l,f)-^+^+^_)P+(l ) t), 

(43) 



dt 

d!L h (n,t) 
dt 



-(p h g +p h c )P h (0,t)+p h d P-(l,t), (44) 



p h c [n h (n-l,t) + P+(n-l,t)} 
(p>*+p'*)Tl h (n,t), forn>2, 



(45) 



^hM. =p h PhM (pj + p fc)n h (i ) t). (46) 

In order to distinguish between the Hill model and our 
Hill-like model with catastrophe-suppressing drugs, we 
replace the subscripts (and superscripts) H of the former 
by h in the latter. 

The steady-state solutions of these kinetic equations 
(see appendix IBl for details) are 



P+(n)=x h (x h + z h ) n - 1 P h (0), 
Ph(n) = y h {x h + z h ) n - 1 P h (Q), 
U^^Zhixh + Zhr-'PhiO), 



where 



with 



and 



Ph(0) 



1 - x h - z h 
l + Vh 



Xh 



p h g {p h d +P h - + )+P h cP- + 
PdiPg + P+-) + PcPd 



Vh 



Pg+Pc 
Pd 



(47) 
(48) 
(49) 

(50) 
(51) 
(52) 



Thus, the equations (HJ)- (J3J) are generalized to the forms 



Zh 



Pc + Pa 



(53) 



dP + (n t) 

— h j t ' = Pgi p h( n -!>*) + n h{n - 1, t)] + pl + P^(n, t) The distributions ([5111531) are decreasing functions of n 

provided (xh + Zh) < 1, i.e, 



(p£+P2+J»+-W(M), forn>2, 



(41) 



(Pg 



Pc) 2 p h - + <p h gP h dP h + _ 



(54) 
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— p„ 


= 0.00 




l>, 


= 0.03 




h 

"" ?c 


= 0.05 




h 

~~ 1\ 


= 0.07 





(a) 



The parameters Pg,Pd, v\- etc. have dimensions 
[time^ 1 ]. For a typical set of values of these param- 
eters, we have plotted the distributions P£(ri), P^{n) 
and Uh(n) in fig[2£a),(b) and (c), respectively, each for 
several different numerical values of pj?. The straight 
lines on the semi-log plots is a reflection of the expo- 
nential decay of the distribution with increasing length 
of the microtubules. Moreover, the longer tails of these 
distributions at at higher values of p^ demonstrates that 
higher p^ causes stronger suppression of the catastrophe 
phenomenon. 




00 



= 1 


; 




— p'= 0.0001 






;/ = 0.001 






— ;i" = 0.01 






~ ;/ = 0.03 






— p* = 0.05 






— ;/ = 0.07 









(c) 



A. Stability of the steady-state 

Let us define the small deviations 

6P±(n)=P±(n,t)-[P±] ss 

SU h (n) = n h (n,t)-[U h ] ss 



(55) 



(56) 



from the corresponding steady-states where the steady- 
state values [P±]ss and [n] ss are given by the equations 
(|17]) , (05]) and P5| , respectively. Linear stability analysis 
of the steady-state distributions leads to the equations 

d5Pt{n,t) = p h [§p + (n _ 1}t) + 5 n h (n - l,t)] 



-(p h q +pt+p*)6P+(n,t) for n>2 (57) 



dSP-(n,t) 



dt 



= p h d 8P^{n + l,t)+p h + _5P+{n,t) 

-(p h d +p-+)sp f ;(n,t) 



for n > 1 (58) 



dSH h (n, t) = p H {mh{n + sp + (n _ 1; t)} 



dt 



{p h g + p^)5n h {n,t)ior n>2 (59) 



FIG. 2: The steady-state distributions (a) P£(n), (b) P^(n) 
and (c) Hh(n) of the size n of microtubules in our Hill-like 
model in the presence of catastrophe-suppressing drugs are 
plotted for several values of p\\ the common parameters are 
Pg = 0.5, pa = 0.75, p+_ = 0.01 and p'l+ = 0.01, each per- 
unit time. 



this condition (|54|) reduces to the condition (fTTj) in the 
limit p h c -> 0. Note that in the limit p h c -> 0, z h -> 
while the expressions (|5ip and ([52")) for Xh and yh re- 
duce to xh and yu given by the expressions §§§ and 
(p~0|) . respectively; hence, in this limit, Hh(n) — > and 
the expressions (jTf)) and (J3HJ) for -P/f( n ) reduce to the 
corresponding expressions © and ([7]) for P^(n). 



d6P+(l,t) _ mh 



dt 



p*SP(0,t)+p1 + SP^(l,t) 
-(p h g +p h c +p h + -)6P+(l,t) (60) 



= p ^P(0,t) (p^+p£)OT*(l,t)> 



Using the normalization condition we get 



(61) 



JP(0, *) = - £ <SP h + (n, t) - £ «V (n, i) - £ ^ (». *) 

n=l 

(62) 



n=l 
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Next, we define 



A± = X>P±(n) 

oo 

A = ^<TL» 



n=l 



and the column vectors 




(63) 
(64) 



Nv 




(66) 



The equations obtained from the linear stability analysis 
above can be written as 



dW(t) 
dt 



M h W(t)+N h <yP^(l,«) (67) 



(65) where the matrix Mh is given by 



-{P H c + P\- - P H a ) 

M h = ( p\_ 




P h -+-P h g 
-P h -+ 

~Pc 



r 



o 
o 

M+Pc) 



(68) 





mi 


ni2 


m 3 





-100.000 


-100.000 


-0.300 


10 


-100.020 


-100.000 


-0.280 


20 


-120.037 


-120.000 


-0.263 


30 


-130.051 


-130.000 


-0.249 


50 


-150.073 


-150.000 


-0.227 


100 


-200.110 


-200.000 


-0.190 


150 


-250.132 


-250.000 


-0.168 



TABLE II: The eigenvalues of the linear stability matrix Mh 
in our Hill-like model. The common parameters are p h g — 
100, p% = 900, P-+ = 0.08, p+_ = 0.22, (each per unit time) 



Note that all the eigenvalues remain negative up to a 
reasonably high value of p c corresponding to the param- 
eter set chosen for the table HU This result indicates that 
the steady state distributions of the lengths of the mi- 
crotubules, which we have derived in this section, remain 
stable up to a moderately high dosage of the catastrophe 
suppressing drug. 



IV. EFFECTS OF DRUGS ON MT: A HYBRID 
HILL-FREED-LIKE MODEL 



In this section we extend the Hill-like model developed 
in section IIIII by taking into account the dependence of 
the rate of growth of the microtubules on the concen- 
tration of the drug-free tubulin subunits in the solution; 
for this purpose we follow the corresponding approach 
developed by Freed (and reviewed in section [TT| for mi- 
crotubule dynamics in the absence of drugs. However, 
the effects of the drug-bound tubulin subunits are taken 
into account in the same way as was done in the Hill-like 
model presented in the section Hill Therefore, the model 



presented in this section is a hybrid of Hill-like and Freed- 
like approaches] the drug-free tubulins are treated fol- 
lowing Freed while the drug-bound tubulins are treated 
following Hill. 

The binding of a drug-bound tubulin subunit with a 
free nucleating site takes place with probability p nc per 
unit time. Following Freed, the concentration of micro- 
tubules of length n which are tipped with drug-bound 
tubulin subunits is denoted by the symbol n(n,i) while 
that of the microtubules tipped with drug-free tubulin 
subunits and in the growing (shrinking) phase is denoted 
by P+(n, t) (P_ (n, t)). Moreover, all the parameters with 
identical susbcripts in this model and in the Freed model 
have the same physical significance. The equations of our 
interest arc 



dP+(n,t) 
dt 



: p gP [P+{n - l,t) + U{n - l,t)} +p-+P-{n, t) 

-(PgP + p c +p + -)P+(n,t) for n >2j39) 



dP_ (n, t) 
~dt 



dU(n,t) 
Jt 



= Pd P-(n + 1A) + p+-P+(n,t) 
(p d +p_+)P_(n,t) for n>l, (70) 



= p c [U(n ~ l,t) + P+(n - l,t)} 
-(PgP + Pc)R{n,t) for n > 2, 



dP+(l,t) 



dt 



■ PnP N + p_ + P_(l,t) 



-(PgP + Pc+P+- )P+(l,t), 



(71) 



(72) 



da(i,t) 

dt 



PncN- { Pg p + p c )n(l,t), (73) 



8 



and 

,00 oc 

-£ = -^pN-p gP Y^[P+(n,t)+n(n,t)]+p d ^2P-(n,t). 

n—1 7i—l 

(74) 

The steady state equations are 

PnP (No - p+ - P- - n) = - PgP (P+ + n) + Pd p_,(75) 

P+{n) = cP+(n - 1) + dP-(n) + cU(n - 1), (76) 



P_(n) = aP-(n - 1) + 6P+(n - 1), 



II(n) = en(n - 1) + eP+(n - 1), 



P+(l)=dP_(l)+/, 



where 



and 



n(i) = .g, 



a = (Pd+P-+)/pd, 



b = -p+-/p d 



c = PgP/(PgP + Pc+P+-), 



d = P-+/{PgP+Pc+P+-), 



e=Pc/(PgP + Pc), 



f = p n pN/(p g p + p c +p + -), 



g =PncN/(p g p + p c ). 



(77) 
(78) 
(79) 
(80) 

(81) 
(82) 
(83) 
(84) 
(85) 
(86) 

(87) 



Note that in the limit of vanishing concentration of drug- 
bound tubulin, i.e., p c — > 0, the expressions (|8ip ~ (j8"4")l 
for a, b, c, and d reduce to the expressions (|26 |) - (f29)) for 
a',b',c' and d', respectively. Moreover, in the limit 
p c — ► and p nc — > 0, equations ([85]) and (|87|) imply 
e — > and g — » while the expression (|86|) reduces to 
the corresponding expression (|30|1 of the original Freed 
model. 
Defining 

00 



^P+(n), 

71=1 



P- =^P-(n), 

71=1 



n = £n(r 



we obtain the following three equations from (|77| . (176)) 
and f75|l . 



(l-a)P_-6P+ = P_(l), 
(1 - c)P+ - cU - dP_ = /, 
(l-e)n-eP+ = n(l), 



(89) 
(90) 



)1) 



We solve for P + , P_ and n using 1(50 ]) , ([9"T] ) and (|75 ) . 
Substituting these solutions in ([89|l . we obtain P-(l) at 
the steady state. Finally, steady state expressions for 
P+(l) and n(l) are obtained from (j?9"j) and (|5U|t . 

The solutions of linear equations (|90|) , (|9"l"j) and ([75]) 
can be obtained in a straightforward way. The solutions 
are 



P+ 



P_ = 



pN PgP 2 (pPn +Pnc) 
V 

pN Pg(pp n +p nc )(pi + PPg) 

V 



n : 



P_(l) = 



where 



N [pp c p n p2 - pp g p- + p nc +PdPncPl] 

V 



(92) 



(93) 



(94) 



pNpPgjpPn + Pnc)iP+-Pd ~ PcP-+ ~ PP-+Pg) 
(PdV) 



(95) 



V = P 3 p 2 gPn + PdPncPl + p 2 p 2 g{Pnc ~ P-+) 
+p 2 PgPn(pi +P2) + PPcP2Pn + PPgPcipnc ~ P- + ) 



with 



and 



+PPgPncP+- + PPgPd(Pnc + P+-), (96) 

Pi=Pc+P+- (97) 



P2 =p d +P-+- 



(98) 



In the following the generating function method is used 
to obtain the P+(n), P_ (n) and n(n) for arbitrary n. We 
proceed by multiplying both sides of equations ((TBI . G7J 
and ([78]) by x n and then suming over n from n = 1 to 
oo. Defining 



and 



(88) 



P+(.T) = ^P + (n) 3 



P_(z) = ]Tp_(n)s 



n(x) = ^n(n)or", 



(99) 



(100) 



(101) 
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we have the following equations for P+(x), P-(x) and 
U(x): 

(1 - cx)P+ (x) - dP- (x) - cxli(x) = xf (102) 



bxP+(x) + (ax - l)P-(x) = -xfi (103) 



exP + (x) + (ex — l)II(a;) = — xg 



(104) 



where (3 = P-(l)- Solving this set of linear equations for 
P+(x), P-(x) and II(x), we find 

P + (x) = 

(aef - acg)x 3 - (fide - eg + ef + af)x 2 + (fid + f)x 



(105) 



P_(x) = 

(beg - bef)x 3 + (bf -fie- fie)x 2 + f3x 



(106) 



U(x) = 

(acg — aef)x 3 + ((3de — eg + ef — ag — bdg)x 2 + gx 
A ■ 
(107) 

where 

A = (ac + ae + bde)x 2 - (a + c + e + bd)x + 1. (108) 



P + (x) = 

(aef - acg)x z - ((3de - eg + ef + af )x 2 + (fid + f )x 



V 2 - 2\y + 1 



P_(x) = 

(beg - bef)x 3 + (bf -/3c- fie)x 2 + fix 
V 2 - 2Ay + 1 



(109) 



U(x) = 

(acg — aef)x 3 + (fide — eg + ef — ag — bdg)x 2 + gx 



y 2 -2Xy + l 



where 



A = (a + e + c + bd)/[2(ac + ae + bde) 1 ' 2 ]. (112) 



and 



y= (a + e + c + bd) 1/2 x. 



(113) 



Solutions for P±(n) and U(n) can be obtained by ex- 
panding term 

(114) 



y 2 -2\y + l 



in the above equations using Taylor series expansion for 
A > 1 and y < 1 and then equating the coefficients of y n 
on both sides. Thus, the expressions for P±(n) and U(n) 
are as follows: 

P+(n) = a^-V/^asUn-^X) - a 2 U n - 2 (\) + aiC/„_ 3 (A)] 

for n > 2, 

(115) 



where 



P+(l) = a 3 U (X) 



a = (ac + ae + bde), 
cti = (aef - acg) /a, 
a 2 = (fide — eg + ef + af)jcfi^ 2 \ 
a 3 = (pd + f), 



and 



V2] 



u n (\) = j2(-iy 



m=0 



(n — m)\ 
m\(n — 2m)! 



(2A) 



n—2m 



(116) 



(117) 
(118) 
(119) 
(120) 



(121) 



with U-i(X) = 0, where the symbol [n/2] represents the 
largest integer smaller than or equal to n/2. Similarly, 

P-(n) = a(" _1 )/ 2 [/J{7„_i(A) + f3 2 U n - 2 (X) + f3iU n - 3 (X)} 

for n > 2, 

(122) 



where 



P_(l) = f3U (X) 



fix = (beg - bef))/a 
fi 2 = (bf -fie-fie)/^/ 2 . 



(123) 

(124) 
(125) 



(110) Finally, 



U(n) - a ("- 1 )/ 2 [ 5 C/„_ 1 (A) + 7 2^- 2 (A) +7i^n-3(A)], 

(126) 



n(l) = gU (X) (127) 



(111) where 



7i = (acg - aef) I a, (128) 
72 = (fide - eg + ef - ag - bdg)/a 1/2 . (129) 

It can be easily checked that these solutions approach 
Freed's solutions in the limit p nc — > 0, p c — > 0. In or- 
der that the steady state solutions for P + (n) , P_ (n) and 
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FIG. 3: The steady-state length distributions (a) P+(n), (b) 
P-(ri) and (c) II(n) of microtubules in our hybrid model for 
several values of p c and compared with corresponding distri- 
butions in our Hill-like model (with appropriate rescaling of 
the parameters and variables following <|13ip ; the common 
parameters are p g = p n — 125, pd = 900, p — \- = 0.08, 
P-t — = 0.22, p nc — pc, and p = 0.8, No = 0.2 (in respec- 
tive units). 



tion (|33|) of the Freed model in the limit p c — > 0. 

In order to plot the steady-state distributions (|1 15[) . 
(|122p and (jl26[) in our hybrid Hill-Freed-like model and 
to compare these with the distributions (|47[) . (|4~5|) and 
for the same set of parameters, we first converted 
the concentrations of the different types of microtubules 
into probabilities and, then, chose the numerical values 
of the parameters so as to satisfy the following relations: 



Pn 



P,I 



P'- + 



= Pd, 
= Pc, 

P+-, 
-P-+- 



(131) 



The distributions plotted in fig[3] shows that with proper 
rescaling of the parameters, as mentioned in (|13ip . the re- 
sults for the Hill-like model and the hybrid Hill-Freed-like 
model are almost identical. Moreover, the higher is the 
value of Pc, the longer are the tails of the distributions; 
this, as explained already in the context of the Hill-like 
model, is a consequence of the stronger suppression of 
the catastrophes by higher p c . 

The mean lengths of the microtubules, which corre- 
spond to the distributions plotted in the fig[3l are plotted 
against p c in figlU Surprisingly, for this set of parame- 
ter values the mean lengths of the microtubules tipped 
with drug-bound tubulin and those tipped with drug-free 
tubulins (both in the growing and shrinking phases) are 
identical for almost all values of p c . However, even for 
this set of parameters values, the fraction of microtubules 
with drug-bound cap increases with increasing p c while 
that with drug- free cap decreases (see figEJ) . 

We now define the "effective" catastrophe frequency 
and the "effective" rescue frequency by the relations 



and 



e// 

P+- 



eff 
P- + 



P+-P+ 



P++P- 



n 



p-+ 



P+ + P- 



n 



P+-P+ 

l-Pn 



P-+P- 

l-Po 



(132) 



(133) 



These "effective" frequencies of catastrophe and rescue 
are plotted against p c in figEl This trend of variation 
is qualitatively similar to the corresponding results of 
laboratory experiments performed with the catastrophe- 
supprssing drug vinblastine (3lj . 



n(n) are decreasing, rather than increasing, functions of 
n, we demand that a < 1 which imposes the following 
constraint on the magnitudes of the parameters: 

P- + {Pc + PVgf < PPgPdP+- (130) 

The condition (|130p reduces to the corresponding condi- 



A. Stability of the steady-state 

We have obtained the exact analytical expressions for 
the matrices that decide the stability of the steady-state 
solutions, which we derived earlier in this section, against 
small deviations. For this purpose we expand the nonlin- 
ear kinetic equations about the steady states and retain 
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FIG. 4: Mean length of the microtubules in the hybrid model. 
The numerical values of all the parameters, which are not 
shown explicitly, are identical to those in figO 



only upto the terms linear in the deviations and drop all 
the terms containing higher orders of deviations. 
Let us define the small deviations 

6p = p(t)-[p] as , (134) 
6P±(n)=P ± (n,t)-[P±} ss (135) 

<m(n)=n(n,i)-[n] ss (136) 

from the corresponding steady-states where the steady- 
state values [P±] ss , [n] ss , etc. are given by the equations 
(fTT5]l . (fTT6)l . Ip2j) . (fl%l) . (fT26l) . (fT27)l . etc. We shall 
also use the symbols [P|_] ss , pP-] S s and to denote 

the steady-state values of P + , P_ and II given by the 
equations 

Expanding the equations (j74"|) about the steady 

state and retaining upto the terms linear in the small 
deviations we get 



dSP+(n) 
dt 



■■ Pg Sp{[P+(n - l)] ss + p(n - 1)] M } 



+P 3 [p] ss {<5P + (« - 1) + <5H(n - 1)} +p_ + 5P_{n) 
-p g Sp[P + (ri)] ss -p g [p] ss SP+(n) 
-(p c +p + -)6P + (n) for n > 2(137) 



FIG. 5: Fractions of the microtubules that are tipped with 
drug-free tubulin and are in the growing phase (P+) or in 
the shrinking phase (P-) and those of microtubules tipped 
with durg-bound tubulin (II). The numerical values of all the 
parameters, which are not shown explicitly, are identical to 
those in the fig[3] 



dSP_ (n) 
dt 



--p d 5P-(n + l)+p + -6P+(ri) 

-{p d +P-+)5P-{n) 

for n > 1 (138) 

Pc {SU(n - 1) + SP+(n - 1)} 



dSn(n) 
Jt 

-Pg[p] ss SIl(n) -p g Sp[IL(n)] sa - Pc SIl(n) 

for n > 2 (139) 



dSP+(l) 



p„T+p_+5P_(l) 



dSH(l) 
dt~ 



-p g {[p] ss SP + (l)+Sp[P + (l)] ss } 

-(Pc+P+-)SP+(l) (140) 

oo oo oo 

= Pnc{- £ SP + (n) 6P -W E ^C")} 

n— 1 n— 1 n— 1 

- Pg {\p]„m(l) + 6p[U{l)] ss } - p c 6U(l) 

(141) 



12 




FIG. 6: Effective catastrophe frequency (a) and effective res- 
cue frequency (b), defined through equations (|132|l and (|I33|) . 
respectively, are plotted against p c for p — 2.0. Note the dif- 
ferent parameters used for the insets. The numerical values of 
all the parameters, which are not shown explicitly, are iden- 
tical to those in the figO 



dS °° °° 

-jjT = -PnT - p g 5p{Y^[P+{n)] ss + ^[n(n)] ss } 



N = 




(147) 



The equations obtained from the linear stability analysis 
above can be written as 



dV(t) 
dt 



MV(i) + N(5P_(1) 



(148) 



n=l 



~Pd 5P_(n)(142) where the matrix M is given by 

n=l 

I 



where 



M = 



P+- 

Pc Pnc 
{Pn -Pg)[p]s 



P-+ -Pn[p]ss 


(Pg -Pn)[p]ss 


M 14 


-P-+ 








Pnc 


-Pg[p]ss ~ Pnc 


-Pg[^]ss 


Pn[p]ss +Pd 


(Pn -Pg)\p]sa 


M 44 


1 

M 44 


= (Pn~Pg)([P+]ss 


+[*]«)- 



(149) 



The characteristic equation is quartic in the eigenvalue 
My = p s [$] ss +p n {N — [P+] ss — [P-]ss — [®]ss} (150) A. We have systematically investigated the trend of vari- 
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p = 0.8M 


Pc 


mi 


ni2 


m 3 


mi 





-99.998 


-0.295 


-22.506 


-100.000 


5 


-105.000 


-0.185 


-22.633 


-104.982 


10 


-110.000 


-0.094 


-22.734 


-109.972 


15 


-115.000 


-0.019 


-22.814 


-114.967 


16 


-116.000 


-0.005 


-22.828 


-115.967 


17 


-117.000 


0.008 


-22.842 


-116.966 


p = 2.0M 


Pc 


mi 


ni2 


m 3 


mi 





-249.999 


-0.299 


-19.566 


-250.000 


50 


-300.000 


-0.121 


-19.717 


-300.023 


100 


-350.000 


-0.030 


-19.785 


-350.051 


126 


-376.000 


-0.001 


-19.803 


-376.061 


127 


-377.000 


0.002 


-19.804 


-377.062 



TABLE III: The eigenvalues of the linear stability matrix M. 

The common parameters are p g — p„ = 125, pd = 900, p |- = 

0.08, p+- = 0.22, p nc = p c , No = 0.2 (in respective units). 



ation of the eigenvalues (7711,7712,7713,7714) of the matrix 
M with the variation of p c for several sets of values 
of the other parameters; those for two values of p are 
given in table iHTl In both these cases, there is a range 
< Pc < p™ a:r of values of p c where all the eigenvalues 
remain negative indicating stability of the steady state. 
Moreover, the higher is the steady-state density p of the 
free tubulins in the solution larger is the value of p™ ax . 
As p c is incresed beyond p" lax , one of the negative eigen- 
values simply changes sign which indicates instability of 
the steady-state. The physical reason for the instabil- 
ity of the steady state at sufficiently high p c is that the 
higher is the p c the larger is the fraction of microtubules 
tipped with drug-bound tubulin which can only grow but 
cannot shrink. 



Thus, although the steady-state distributions of the 
microtubules in the Hill-like model and those in the hy- 
brid Hill- Freed- like model can be made practically identi- 
cal by making appropriate correspondence of the param- 
eters in the two models, the latter has a richer dynamics 
which is reflected also in the linear stability analysis. 



V. SUMMARY AND CONCLUSION 

In this paper we have developed models for study- 
ing some generic effects of a class of drugs on 
the polymerization-depolymerization dynamics of micro- 
tubules in the absence of GTP and GDP. The class of 
generic drugs under consideration suppress catastrophes; 
more specifically, we assumed that (i) the microtubules 
capped with drug-bound tubulins do not exhibit catas- 
trophe, and (ii) the rate constant for the attachment 



of a drug-bound tubulin is, in general, different from 
that of a drug-free tubulin. Although the effects of real 
catastrophe-suppressing drugs are known [26[ to be more 
complicated than those of the generic model assumed 
here, some predictions of the model are in good quali- 
tative agreement with the corresponding effects of vin- 
blastine over a limited parameter regime. 

One of the two models, namely, the Hill-like model 
proposed here is an extension of the model devel- 
oped by Hill [| to describe microtubule polymerization- 
depolymerization dynamics in the absence of drugs. Al- 
though mathematical treatment of this model is quite 
simple, it does not take explicitly account for the dynam- 
ics of the tubulin concentration in the solution. There- 
fore, we have also developed a more detailed model in 
which the effects of the concentration of the pure (i.e., 
drug-free) tubulin subunits are taken into account in our 
model in a manner similar to that done by Freed [2(jj in 
his recent theoretical study of the microtubule dynamics 
in the absence of drugs. However, in this model the ef- 
fects of the tubulin subunits bound to the drug molecules 
could not be taken into account in a similar manner. 

The reason for the difficulty in treating the concen- 
trations of the drug-free and drug-bound tubulins in so- 
lution on an equal footing is generic to the Hill-Freed 
approach. In order that the dynamical equation for the 
concentration of drug-bound tubulin subunits can reach 
a steady-state, we have to allow both attachment as well 
as detachment of the drug-bound subunits to the micro- 
tubule. However, if shrinking of a microtubule tipped 
with drug-bound tubulin is allowed, then, immediately 
after such a detachment, one needs to know the status of 
the subunit at the newly exposed tip, i.e., whether the 
tip consists of a drug-free tubulin or a drug-bound tubu- 
lin. But, in the Hill-Freed approach, the system does 
not have a memory of the past history in the sense that 
the the model does not keep a record of the status (i.e., 
whether or not bound to a drug molecule) of all the tubu- 
lin subunits, starting from the nucleation center up to the 
tip. 

Therefore, in order to overcome this technical diffi- 
culty, we have incorporated the effects of the drug-bound 
tubulin subunits in a manner similar to the approach fol- 
lowed in the Hill model 9]. Thus, our second model may 
be regarded as a hybrid of the Hill-like and Freed-like 
modeling strategies for the concentrations of the tubulin 
subunits. 

For both the Hill-like model and the Hill-Freed-like 
hybrid model we have derived exact analytical expres- 
sions for the steady-state probability distributions of the 
lengths of microtubules tipped with drug-bound tubu- 
lin subunits as well as those of microtubules tipped with 
pure (i.e., drug-free) tubulin subunits in the growing and 
shrinking phases. We have also compared the trends of 
variations of some of the relevant quantities with the vari- 
ation of the dosage of the catastrophe-suppressing drug. 

We have carried out linear stability analysis of the 
steady-states and established that in both the models 
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the length distributions of the microtubules remain sta- 
ble unless p c becomes sufficiently high to destabilize the 
steady-state. 



where x and z are given by the equations (|5Tj) and (j53[) . 
respectively. The special case corresponding to n = 1 of 
the equation dPtn,t)/dt = 0, leads to the steady state 
solution for P_(2), 



APPENDIX A: ON THE USE OF CHEBYSHEV 
POLYNOMIAL 



It is straight forward to see that 
a' + c' + b'd' = 1 



PgPiPd +P F +) 



Pd(PgP + P+) 



= 1 + a'c' 
Therefore, re-expressing A' as 



1 + old 



1 



2(a'c') 1 / 2 2 
we find, in general, A' > 1 



|(aV) 1 / 4 -(aV)- 1 / 4 | 2 + 2 



(Al) 



(A2) 



APPENDIX B: STEADY-STATE SOLUTIONS OF 
HILL-LIKE MODEL WITH DRUGS 



P_(2) = —[(p d + p-+)y - p+-x]P(0). (B4) 

Pd 

(B5) 



A little bit of algebraic manipulqation leads to 



P_(2) =y(z + z)P(0). 



(B6) 



Steady-state solutions for P+(2) and 11(2) can be ob- 
tained in a similar way. Steady-state solutions for P_ (n), 
P+(n) and n(n), for arbitrary n, are straightforward gen- 
eralizations of our observations upto n = 4. 

Since P+(n), P-(n) and H(n) are probabilities of find- 
ing microtubules of n subunits in growing, shrinking and 
in catastrophe-arrested phase, we expect the following 
normalization condition 



Here we obtain the steady-state solutions of the equa- 
tions (|4"Tj) - (|4T))) for Hill- like model in the presence of drugs 
by equating the RHS of these equations to zero. The solu- 
tion for P-(l) is obtained by demanding dP(Q, t)/dt = 0. 
This leads to 



P_(l) = 



(Pg 



Pc) 



Pd 



P(0) = yP(0) 



(Bl) 



where y is given by the equation (|52|) . Similarly, claiming 
dP + (l,t)/dt = and «fll(l, t)/dt = 0, we have 



H(l): 



[PgPd +p- + {p g +Pc)\ 



Pd{Pg +P, 
P. 



Pc+Pg 



P+-) 

P(0) = zP(0). 



P(0) = xP(0)(B2) 
(B3) 



P(0) + J2(P+(n) + P-(n) + n(n)) = 1. (B7) 



This leads to the form ([HI]) for P(0). 
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